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X. 

We analyze the spectral moments of the V-A two-point correlation function. Us- 
ing all known short-distance constraints and the most recent experimental data from 
tau decays, we determine the lowest spectral moments, trying to assess the uncer- 
tainties associated with the so-called violations of quark-hadron duality. We have 
generated a large number of "acceptable" spectral functions, satisfying all conditions, 
and have used them to extract the wanted hadronic parameters through a careful sta- 
tistical analysis. We obtain accurate values for the xPT couplings L\q and C$7, and 
a realistic determination of the dimension six and eight contributions in the operator 
product expansion, C 6 = (-5.4 + • 10~ 3 GeV 6 and O s = (-8.9 ± \ 2 f) ■ 10~ 3 GeV 8 , 
showing that the duality-violation effects have been underestimated in previous lit- 
erature. 



1 Introduction 



QCD sum rules (QCDSRs) [HE] have been widely used during the last thirty years to study 
many important aspects of QCD. They constitute a very useful tool, enabling us with a 
powerful connection between QCD parameters and physical observables. 

The basic assumption behind the QCDSR techniques is that the quark and hadron 
degrees of freedom provide two dual descriptions of the same strong interaction dynamics. 
This quark-hadron duality is a consequence of the assumed confinement of QCD. In more 
technical terms, a QCDSR is a dispersion relation relating the value of a given two-point 
correlation function at some Euclidean value of Q 2 with an integral over the corresponding 
spectral function in the Minkowskian domain. Quark-hadron duality allows us to calculate 
this Minkowskian integral in terms of hadrons, using the available experimental data. 
Ideally, the resulting QCDSR is an exact mathematical relation arising from analyticity and 
confinement (duality). In practice, however, a series of approximations need unavoidably to 
be adopted in its specific numerical implementation. In the Euclidean region the correlator 
is approximated by its short- distance Operator Product Expansion (OPE) [3], truncated to 
a given finite order. In the Minkowskian region, since experimental data is only available 
at low energies, the integral over the physical spectral function is usually cut at a certain 
finite invariant-mass Sq; from so up to oo, one then adopts the short- distance information 
provided by the OPE. 

The uncertainties associated with all these approximations are usually known as vio- 
lations of quark-hadron duality. They are difficult to estimate, because of our inability 
to make reliable QCD calculations at low and intermediate energies. The normal way to 
assess the theoretical uncertainties of QCDSRs consists in estimating the OPE truncation 
error and testing the stability of the results with variations of so- However, this method 
is too naive and can underestimate the effects not included in the OPE, i.e. the difference 
between the physical correlator and its OPE approximation. 

Violations of QCD quark-hadron duality [I] have been relatively poorly studied and 
often disregarded. Its importance in finite energy sum rules (FESRs) has attracted some 
attention recently [SHE], owing to the phenomenological need for higher accuracies. To 
estimate the size of these effects is of course of maximal importance, if we want to master 
the strong interaction at all energies and be able to perform precision QCD calculations. 
This importance extends to all particle physics when one realizes that those calculations are 
often necessary to disentangle new physics from the Standard Model. Moreover, duality vi- 
olations will also be present in new-physics scenarios characterized by a strongly-interacting 
dynamics. A better knowledge of duality violations in QCD would help to understand their 
role in more exotic theories. 

In the following, we present a detailed analysis of the possible numerical impact of 
duality violations in the description of the two-point correlation function of a left- and a 
right-handed vector currents. This is a very good laboratory to test the problem because 
this correlator is an order parameter of chiral-symmetry breaking: in the massless quark 
limit it vanishes to all orders in perturbation theory; its operator product expansion only 
contains power-suppressed contributions, starting with dimension six. In the absence of any 
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Figure 1: Analytic structure of II(s). 



theory of duality violations, we will use a generic, but theoretically motivated, model [HE] 
to assess the phenomenological relevance of these effects. 

The theoretical ingredients of our analysis are presented in the next section. Section [3] 
contains a detailed discussion of the behaviour of the physical spectral function at high 
energies. Using the most recent experimental data, we generate a large number of "accept- 
able" spectral functions which satisfy all known QCD constraints. Our numerical results, 
obtained through a careful statistical analysis of the whole set of possible spectral functions, 
are given in section |H Section [5] summarizes our findings. 

2 Theoretical Framework 

The basic objects of the theoretical analysis are the two-point correlation functions of the 
vector and axial- vector quark currents J^[x\ defined as follows: 

n£» = ^|d 4 xe^(0|T(^(x)^(0)t)|0) 

= + ) n ^(? 2 ) + <?<t n^ 2 ) • (i) 

Although our analysis can be applied to any correlation function, we will only study here 
the non-strange correlators and therefore J^ix) will denote the Cabibbo-allowed vector or 
axial- vector currents, V^ d (x) = uj^d and A^ d = uj^^d. Moreover, we will concentrate on 
the J = 0+1 part of the V — A difference, that is nothing but the correlation function of the 
left- and right-handed currents, L» d (x) = V^{x) - A» d (x) and R^Jx) = V£ d (x) + A» d (x), 
that is 

n( s ) = n^\s)-u^(s) = -M^ + n( s ), (2) 

where we have made explicit the contribution of the pion pole to the longitudinal axial- 
vector two-point function. We will work in the isospin limit, m u = m^, where y(g 2 ) = 0. 
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The correlator Tl(s) is analytic in the entire complex s-plane, except for a cut on the 
positive real axis which starts at the threshold s t h — 4m^. Applying Cauchy's theorem in 
the circuit in Fig. [T]to the function w(s) n(s), one gets the exact relation: 

/ ds w(s) p(s) + / dsw(s)U(s) = 2f%w(ml)+Res[w(s)Il(s)], (3) 
Js th 2 ™ J lsl=SQ ■ " s=o 

where p(s) = -Imll(s), w(s) is a general analytic weight function except maybe at the 
origin where it can have poles, and Res [-F(s)] is the residue of F(s) at s — 0. Integrals of 

s=0 

the chiral spectral function p(s) times w(s) from threshold s t h up to s are usually called 
spectral chiral moments Mw s )(so); when so — >■ oo we will denote them M w ( s ) for brevity. 

In order to evaluate the contour integral of ([3]), one approximates H(s) with its OPE 
expression 

n op E(s) = ^ www s ^ ft* , (4) 

where (02k){v) are vacuum expectation values of operators with dimension d = 2fc; their 
associated Wilson coefficients C 2 kiy) contain logarithmic dependences with — s. Notice 
that both are //-dependent quantities, but this dependence cancels in their product 2 k- 
The use of the OPE introduces a systematic error in the relation ()3]), which is expected to 
be dominated by the region close to the positive real axis where the OPE approximation 
does not apply (except at sq = oo [10J). 



Let us rewrite Eq. ((3]) in the form: 

f°dsw(s)p(s) + — I dsw(s)n OPE ( 
Js th 2vn J ls{=so 



s) + DV[w(s),s ] 



2/>K) + Res[ W ( S )n( S )], (5) 

s=0 



where 



DV[w(s), s ] = -L / ds w(s) (n( a ) - n OPE ( s )) 



(6) 



parameterizes the violation of quark-hadron duality that we are interested in. Notice 
that DV[w(s),s ] depends on the weight function w(s) and on the circuit-radius s . The 
relation (J5]) contains all the elements of a standard sum rule. The first term is the hadronic 
part, that in our case is nothing but an integral of the V—A non-strange spectral function 
that has been measured in r decays (for s < ml) [m - [T5] . while the second term is the OPE 
contribution to the contour integral at \s\ = sq. The second line contains the pion-pole 
contribution and the residue at the origin for negative power weight functions, l/s n , which 
is calculable with Chiral Perturbation Theory (%PT) [IB] . 

Sum rules of this type have been applied countless times in the last thirty years in 
order to extract theoretical parameters like quark masses [I7J[T8], the strong coupling 
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constant [19], QCD condensates [SlUMU or ^PT couplings [IB [23 [IMS]. Or, used in 
the other way around, to make predictions of hadronic observables. 

In the chiral limit (m u = = 0) the correlator Tl(s) vanishes identically to all orders 
in perturbation theory and therefore its OPE contains only power-suppressed contributions 
from dimension d = 2k operators, starting at d — 6, as we have already indicated in (pO). 
The nonzero up and down quark masses induce tiny corrections with dimensions two and 
four, which are negligible at high values of s. This makes this correlator a very interesting 
object in the study of non-perturbative QCD. 

In order to analyse duality violation (DV) effects in different sum rules, we will use the 
weights w(s) = s n , with n = —2, —1, 2, 3, that generate the following four FESRsS 

M_ 2 (s ) = / ds-p(s) = 16C 8 f - DV[l/s 2 ,s ], (7) 
f s ° 1 

M_i(s„) = / ds-p(s) = -81$ - DV[l/s,s ], (8) 

/•so 

M 2 (s ) = / dss 2 p(s) = 2/>* + O, - DV[s 2 ,s ], (9) 
J «th 

/•SO 

M 3 (s ) = / dss 3 p(s) = 2/X - °s ~ DV[s 3 ,s ], (10) 

Js th 

where L^q = — |n(0) and Cgf = ^ II (0) are quantities that can be written in terms of 
low-energy constants [HI], while Oq^ are defined in Eq. (@J. These four sum rules 
have been used in the past [T3 1I20H32] to extract the values of either the couplings 
Lio and Cs7, or the vacuum expectation values of the dimension six and eight operators 
appearing in the OPE. In those works the DV effects were just inferred from the so-stability 
(if not just neglected), that as we will see can be a misleading method. Here we want to 
analyze the effect of DV on these four observables using a different approach that will be 
introduced in the following sections. 

For the computation of the hadronic integral representation of the moments M n (so) 
we will use the 2005 ALEPH data on semileptonic r decays [H], shown in Fig. |2l which 
provide the most recent and precise measurement of the V — A spectral function p(s). 



2.1 Theoretically-known spectral moments 

In the four sum rules introduced in the previous section, we use the experimental data 
to extract theoretical information, namely the value of the corresponding parameters or, 
equivalently, the value of the spectral moments for sq — > oo, M n . There exist a few 
additional sum rules where we know theoretically the value of the spectral moments when 
So — > oo. These sum rules will play a special role in our analysis because they give us very 

1 Here we neglect the logarithmic corrections to the Wilson coefficients in the OPE. The error associated 
to this approximation is expected to be smaller than the other errors involved in the analysis, as was found 
e.g. in Refs. [lfl[33]. 
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Figure 2: Non-strange V-A spectral function p(s) = -lmU^y_ A (s) measured from 
hadronic r decays by ALEPH [TT] . 

valuable information on the spectral function p(s) for s > Sq. The three sum rules that we 
will use are: 

/•OO 

M = / dsp(s) = 2fl, (11) 



oo 



/oo 
dssp(s) = 2/X, (12) 
-th 

/ ds s log (-4) p(s)L g=0 = W-m^EM^/o 2 . (13) 

The relations f llip and (1121) are the well-known first and second Weinberg sum rules 
(WSRs), while the third identity is the pion sum rule (7rSR) giving the electromagnetic 
pion mass splitting in the chiral limit [35]. In the second WSR there are contributions of 
the form O{m 2 q asso) |36j, where so is the upper limit of the integral, but they are negligible 
for the values of Sq that we are considering. 



2.2 Duality violation 

To get vanishing DV in sum rules like fl5]) and (EHTUJ) one could think working with an 
infinite Cauchy radius so, but this is clearly not an option because the spectral function p(s) 
is only known up to s max = m 2 .. We can predict the value of p(s) at high-enough energies 
using perturbative QCD, but there is an intermediate region above s max where perturbation 
theory is still not reliable. Therefore we have to deal with this DV unavoidably, and it is 
important to keep in mind that at so ~ 3 GeV 2 it can represent a sizable contribution to 
the sum rules, as the WSRs show clearly (see e.g. Fig. 1 in Ref. [22J). 

Since the solution to QCD is not known yet, DV is almost by definition a non-calculable 
quantity and that is the reason why it has been taken to be negligible very often. But in 
order to make precise and reliable predictions one must worry about the size of this effect. 
As it is commonly done, we have defined the DV in Eq. as the uncertainty associated 
with the use of the OPE. The usual strategy to estimate the size of the DV has been to 
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look at the stability of the spectral moments with variations of sq. This stability can be 
improved adopting the so-called "pinched weights" [23], polynomial weight functions with 
a zero at s = so which suppresses the contribution to the integral (El) from the the region 
close to the positive real axis. As we will see, the stability with so obtained with these 
weights can be misleading in some situations. 



Taking into account that DV[w(s), s ] vanishes for s — > oo, one can easily re-write 
Eq. © in the form [SHZIE] 



expressing the DV effect as an hadronic integral that can be analyzed phenomenologically. 

We know from QCD that the spectral function p(s) has to vanish at high values of s 
and, consequently, we expect the region right above s to be the most relevant in ( JT4")) . 
This makes the "pinched weights" an interesting tool to minimize the DV. However, in 
(TH1) we can see something that is hidden in ((6]), namely that one has to worry also about 
the possible enhancement of the contribution from the high-energy part of the integral 
(s so) produced by the "pinched weights". And thus, we see that the use of these 
weights can worsen the situation. Another direct consequence from (THl) . unless accidental 
cancelations occur, is that by weighting less the high-energy part of the spectral integral 
one can get smaller DV. In particular, for our spectral moments M n (so), one expects the 
DV effects to increase with increasing values of n. Thus, the size of the DV will be smaller 
in the determination of L\q than in the determination of the chiral moment M 2 . 

To quantify the DV uncertainties of a given sum rule we must then estimate the possible 
behavior of the spectral function beyond so- The DV is an estimate of the freedom in the 
behavior of the spectral function above so, once all the theoretical and phenomenological 
knowledge on that spectral function and on its moments has been taken into account. For 
instance, QCD tells us that p(s) must go quickly enough to zero when s — > oo. This is 
a valuable information, but one can still imagine infinite possible shapes for the spectral 
function and, therefore, the limits imposed on DV effects are poor and not good enough 
for most phenomenological analyses. 

Some theoretically motivated models for the DV were advocated in Ref. [3j. We will 
adopt a simple parameterization of the spectral function at high energies, based in the 
resonance model proposed in [I] and similar to the one used in Refs. 0[9]. Following the 
discussion above, we add more physical constraints to the behaviour of p(s) and require 
that it satisfies the WSRs and the 7rSR [6] . Our goal is to generate a bunch of physically 
acceptable spectral functions and translate this information into DV limits. 

A similar work has been done in [TIE] to estimate the DV uncertainties associated with 
the determination of a s from hadronic r decay data. An important difference of our present 
study with those works is that they make separate analyses for the vector and axial-vector 
channels, without imposing the constraints from the WSRs and 7rSR. In fact, one can easily 
check that those sum rules are not satisfied for the vast majority of the generated spectral 
functions used in [312] (as can be seen in Fig. 2 of ref. |Sj). So the results found there 
cannot be applied to the V — A channel that we want to study here. 




(14) 
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3 Acceptable V — A Spectral Functions 



3.1 Spectral-function parameterization 

We split the integral of the spectral function p(s) in two parts. For the low-energy part of 
the integral we will use the ALEPH data, whereas in the rest of the integration range we 
will work under the assumption that the spectral function is well described by the following 
parameterization 

p(s > s z ) = k e" 7S sin(/3(s - s z )) , (15) 

that has k, 7, /3 and s z as free parameters. From the ALEPH data we know that the V — A 
spectral function p(s) has a second zero around 2 GeV 2 (see Fig. EJ), which is represented 
in our parameterization through the s z parameter. We will take this zero as the separation 
point between the use of the data and the use of the model. 

At high values of s this parameterization appears naturally in the equidistant resonance- 
based model with finite widths introduced in [1] . It has also been used for the vector and 
axial- vector correlators in Ref. j9], based on the expected exponential fall-off associated 
with the intrinsic error of an asymptotic expansion; the sine function reflects the periodicity 
of the daughter trajectories in the spectrum of the Regge theory. 

In the region 2.0 GeV 2 < s < 3.3 GeV 2 the proposed parameterization is compatible 
with the ALEPH data; the corresponding x 2 fit gives the result^] 

xL„(«, 7, A s z ) = x 2 (l-00, 1.05, 0.40, 2.03) = 4.4 « d.o.f. = 43 . (16) 

In fact the compatibility appears to be too good, in the sense that the minimum x 2 is much 
smaller than the number of degrees of freedom (d.o.f.): 43 = 45 points - 2 parameters. This 
low value of Xmin was a l so found in Refs. [91137] . 

3.2 Imposing constraints 

As we have already said, the WSRs and the 7rSR in (ITTj) . (fT2l and fflBl are an important 
source of information on p(s), for s values beyond the range of the r data. In the literature, 
the use of this information has been mostly limited to define the so-called "duality points" , 
values of s for which the WSRs are satisfied, i.e. DV[s n , s°] = (n — 0, 1). These duality 
points are frequently used to evaluate the other FESRs, but this introduces an unknown 
systematic error and several ambiguities, like which duality point is the best option. 

We will fully use that information by imposing that the spectral function p(s), given 
by the latest ALEPH data below s z ~ 2 GeV 2 and Eq. (fl5l) for s > s z , fulfils the two 
WSRs and the 7rSR within uncertainties. This requirement constrains the regions in the 
parameter space of model (fT5l) that are compatible with both QCD and the data. We will 
find all possible tuples^ (k, 7, /3, s z ) which are compatible with such constraints by fitting 

2 Hereafter, unless otherwise stated, we include all correlations among the points. 
3 We will talk about "tuple" referring to a set of values (k,7,/3, s z ). 
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the model. In this way, we analyse how much freedom is left for the shape of the spectral 
function after imposing all we know on p(s) from data plus QCD. We will also require the 
compatibility between models and data in the region^ 1.7 GeV 2 < s < 3.15 GeV 2 . 
The four imposed conditions can be written quantitatively in the following formal 

PS Z POO 

/ p(s) ALEPH ds + / p(a;«, 7 ,/3,s 4 ) ds = 2 ft = (17.1 ± 0.4) ■ 10" 3 GeV 2 , (17) 

J0 J s z 

ps z poo 

/ p(s) ALEPH s ds + / p(s',K,j,P,s z ) s ds = 2/ 2 m 2 = (0.3 ± 0.8) • lO" 3 GeV 4 , (18) 

JO Js z 

/ p(s) ALEPH s log ( — - — ^ ) ds + p(s: k, 7, B. s z ) s log ( — - — ^ ) ds 
Jo FV ' VlGeV 2 / Js z VIGeV 2 / 

= (m 2 - m 2 + ) EM fl = -(10.9 ± 1.5) ■ 10" 3 GeV 4 , (19) 
X 2 (^,6,s z ) < xL = d.o.f. = 54. (20) 



3.3 Selection process of acceptable models 

After defining the minimal conditions that a tuple has to satisfy in order to be accepted, 
we perform a scanning over the 4-dimensional parameter space, looking for physically 
acceptable tuples. We emphasize the importance of taking properly into account the data 
correlations. For instance, if one analyses the compatibility of a null spectral function with 
the ALEPH data in the region (2, 3.15) GeV 2 , the resulting minimum \ 2 is very sensitive 
to these correlations: 

X 2 (0.0,7,/3,s z )/d.o.f. = 0.99 (correlations included), (21) 
X 2 (0.0,7,/3, s 2 )/d.o.f. = 4.58 (correlations excluded) . (22) 

To perform the parameter-space scanning process, we adopt the following procedure. 
First, we define a rectangular region such that it contains the four-dimensional ellipsoid 
defined by x 2 (/t, 7, B, s z ) = d.o.f., and we create a lattice with 20 4 = 16 • 10 4 points, that 
is, 16 • 10 4 tuples (or functions). We find that 1789 of them satisfy our set of minimal 

4 Although we are assuming that the model describes correctly the spectral function beyond s z ~ 2 GeV 2 , 
we impose the compatibility with the data from 1.7 GeV 2 to ensure the continuity of the spectral function 
in the matching region between the data and the model. 

5 The quoted errors in Eqs. (|17l) and (|18[) are just data errors, whereas in (|19p the main uncertainty 
comes from the fact that quark masses do not vanish in nature and we are using real data (not chiral-limit 
data). We estimate this uncertainty taking for the pion decay constant the value fo = 87 ± 5 MeV, that 
covers a range that includes the physical value and the different estimates of the chiral limit value [38] . 
We also include a small uncertainty coming from the residual scale dependence of the logarithm, which 
is proportional to the second WSR. We consider X ~ 1 GeV a good choice of scale because higher values 
would suppress the high-energy part of the integral (the information that we want to use), while smaller 
values would generate larger r-data errors in (|19[) . losing also information about the high-energy region. 
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Figure 3: Statistical distribution of acceptable models in the parameter space k (upper- 
left), 7 (upper-right), (lower-left) and s z (lower-right). 

conditions; i.e., 1789 of them represent possible shapes of the physical spectral function 
beyond 2 GeV 2 . Fig. [3] shows the statistical distribution of the parameters of our model 
after the selection process. In Fig. H]we show the distribution of the quantity x 2 ( K > 7> A s z ) 
for those tuples that have passed the selection process. We find that all accepted tuples 
generate values of \ 2 larger than 10.0; i.e., tuples following the central values of the ex- 
perimental points do not pass the selection process; neither do the tuples that go above 
the central values. Thus our model indicates clearly that the third bump of the spectral 
function should be smaller than what the ALEPH data suggest (see Fig. [2]). The size of 
this third bump is an important issue that future high-quality r decay data could clarify. 
For illustrative purposes, Fig. [5] shows one of the hundreds of functions that satisfy our set 
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Figure 4: Distribution of x 2 ( K ? 7? P, s z ) values for acceptable tuples. 
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of conditions. 




Figure 5: Spectral function p(s) generated with (k,j,{3,s z ) = (0.24,1.23,2.82,2.03), to- 
gether with the experimental ALEPH data [TT]. \ 2 = 38.7 for this tuple. 



4 Numerical Results 

For each one of the hundreds of functions that have passed our selection process, we can 
calculate the associated values of Cfj? , L\q, 0% and Og, simply carrying out the integrals 
of Eqs. (TT HTOj) with s — > oo. The results of this analysis are summarized in Fig. El which 
shows the statistical distribution of the calculated parameters. From these distributions, 
one gets the final numbers: 

C§ = (8.167t£^ ± 0.12) • 10" 3 GeV~ 2 = (8.17 ± 0.12) • 10~ 3 GeV~ 2 , (23) 
Lf = (-6.46 ^;° 3 ± 0.07) -HT 3 = (-6.461°$) • lO" 3 , (24) 
C> 6 = (-5.4 + 3 ;<5 ± 1.2) ■ 10~ 3 GeV 6 = (-5.4 1 3 ; 6 ) ■ 10~ 3 GeV 6 , (25) 

8 = (-8.9^ 2 1 4 ±2.l)-10- 3 GeV 8 = (-8.9 ± \ 2 f) ■ 10~ 3 GeV 8 , (26) 

where the first error is that associated to the high-energy region (integral from s z to infin- 
ity) , that we compute from the dispersion of the histograms of Fig. El and the second error 
is that associated to the low-energy region (integral from zero to s z ), that we compute in a 
standard way from the ALEPH data. This results correspond to the 68% probability region 
(one sigma). Since the first error is not gaussian we show also now the 95% probability 
results (95% of the acceptable spectral functions give a result within the quoted interval): 

c$ = (8.167 ±8:oot ± °- 24 ) • 10-3 GeV ~ 2 = ( 8 - 17 ± °- 24 ) • 10-3 GeV ~ 2 » ( 27 ) 

Lf = (-6.46±^±0.14)-10- 3 = (-6.46l8i8 • 10" 3 , (28) 
C> 6 = (-5.4l| 2 ±2.4) • 10~ 3 GeV 6 = (-5.4 ±£*) ■ 10" 3 GeV 6 , (29) 
O s = (-8.9 1 Hi ± 4.2) • 10~ 3 GeV 8 = (-8.9 1 \U) • 10" 3 GeV 8 . (30) 
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Figure 6: Statistical distribution of values of Cg (upper-left), Uf^ (upper-right), Oq (lower- 
left) and O s (lower-right) for acceptable models. 

Our calculations have been done with a very simple, but physically motivated, param- 
eterization of DV [U|9]. Most likely this parameterization does not represent the actual 
shape of the V— A spectral function, but it accounts for the possible freedom of the function 
p(s) beyond 2 GeV 2 and its consequences on the observables. Our statistical analysis trans- 
lates the present ignorance on the high-energy behaviour of p(s) into a clear quantitative 
assessment on the uncertainties of the phenomenologically extracted parameters. 

As expected, the DV effects have very little impact on the values of Cff and L^, 
because the corresponding FESRs (j7]) and (jSj) are dominated by the low-energy region 
where the available data sits. Our results are in excellent agreement with the most recent 
determination of these parameters, using the same ALEPH r data, performed in Ref. |31j : 
Cg = (8.18 ± 0.14) • 1(T 3 GeV" 2 and Lf Q = -(6.48 ± 0.06) • 10~ 3 . 

The situation is not so good for the moments M 2 and M 3 (or equivalently Oq and O s ), 
which are sensitive to the high-energy behaviour of the spectral function. The present 
ALEPH data, together with the constraints from the WSRs and the 7rSR, are not good 
enough to determine the sign of Og] the DV uncertainties turn out to be too large in this 
case. Our results are slightly better for Oq, where there is no doubt in the sign, but again 
the effects of DV imply larger uncertainties than what was estimated in previous works. 
Our results are compared in Fig. [7] with previous determinations of Oq and Os- One 
recognizes in the figure the existence of two groups of results that disagree between them. 
For Oq there is a small tension between a bigger or smaller value, whereas in the case of 0% 
the disagreement affects to the sign and is more sizeable. In some cases the discrepancy 
appears to be related with a two-fold ambiguity in the adopted choice of "duality points" . 
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Figure 7: Comparison of our results for Oq (left) and Og (right) with previous determina- 
tions IIIHIllriE^lElEZMIMfll] (we show for every method the most recent determi- 
nation). The blue bands show our results at 65% C.L., while the 95% probability regions 
are indicated by the dotted lines. 

Our analysis indicates that the DV error was grossly underestimated in most of the previous 
works based on FESRs ([7]) - ({TO]) . Only Refs. [22j[28] quote uncertainties similar to ours, 
although our error bands are slightly shifted in such a way that the tension with the other 
estimates is reduced. 



5 Summary 

The phenomenological requirement for increasing precisions in the determinations of ha- 
dronic parameters makes necessary to assess the size of small effects which previously 
could be considered negligible. In particular, a substantial improvement of QCDSR results, 
needed to determine many hadronic observables both in the Standard Model and in models 
beyond it, could only be possible with a better control of DV. 

Violations of quark-hadron duality are difficult to estimate because those effects are 
unknown by definition. They originate in the uncertainties associated with the use of the 
OPE to approximate the exact physical correlator. As defined in Eq. ([6]), DV effects corre- 
spond to an OPE approximation performed in the complex plane, outside the Minkowskian 
region, which deteriorates in the vicinity of the real axis. Using analyticity, the size of DV 
can be related with an integral of the hadronic spectral function from sq up to oo, given 
in Eq. ffl^j) . which allows us to perform a phenomenological analysis. 

We have studied the possible role of DV in the two-point correlation function n(s). 
This V — A non-strange correlator is very well suited for this analysis because: i) it is a 
purely non-perturbative quantity in the chiral limit, ii) there are well-known theoretical 
constraints, and iii) there exist good available data from r decays. Moreover, different 
moments of its spectral function provide hadronic parameters of high phenomenological 
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relevance. 

We have assumed a generic, but theoretically motivated, behaviour of the spectral func- 
tion at high energies, where data are not available, with four free parameters. This allows 
us to study how much freedom in p(s) could be tolerated, beyond the requirement that 
all known QCD constraints are satisfied. Performing a numerical scanning over the four- 
dimensional parameter space, we have generated a large number of "acceptable" spectral 
functions, satisfying all conditions, and have used them to extract the wanted hadronic 
parameters through a careful statistical analysis. The dispersion of the numerical results 
provides then a good quantitative assessment of the actual uncertainties. 

We have determined four hadronic parameters of special interest: Cff , L^, Oq and 
Og. Our final numerical results are given in Eqs. (!23|) -( l26|) for the one sigma results and 
( p7|) -( |30|) for the 95 % probability results. The parameters and are in excellent 
agreement with the most recent determination using FESRs and the same ALEPH data 
|31j . The vacuum condensate Oq is an important input for the calculation of the CP- 
violating kaon parameter e' K , it dominates the A J = 3/2 contribution to e' K [2"2~||2"3"]. The 
determination of this contribution is an important goal of lattice QCD calculations and 
independent information is required to test the reliability of those results. We will study 
the consequences of our results for e' K in a forthcoming publication [42J. 

Our analysis indicates that the DV error was grossly underestimated in most of the 
previous QCDSR determinations of Oq and Og based on the FESRs (I7j) - ( TTU1) . The present 
V—A non-strange tau data between 2 GeV 2 and 3 GeV 2 [H] is not good enough to constrain 
the spectral function with the needed accuracy. Good data in that energy region with much 
smaller experimental uncertainties is clearly required. Future high-statistics r-decay data 
samples could allow a substantial improvement of our results, helping to clarify the actual 
high-energy behaviour of the V — A spectral function. 
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